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We derive a coarse grained, free-energy functional which 
describes droplet configurations arising on nucleation of a 
product crystal within a parent. This involves a new 'slow' 
vacancy mode that lives at the parent-product interface. A 
mode-coupling theory suggests that a slow quench from the 
parent phase produces an equilibrium product, while a fast 
quench produces a metastable martensite. In two dimensions, 
the martensite nuclei grow as 'lens-shaped' strips having alter- 
nating twin domains, with well-defined front velocities. Sev- 
eral empirically known structural and kinetic relations drop 
out naturally from our theory. 

PACS: 81.30.Kf, 81.30.-t, 64.70.Kb, 64.60.Qb, 63.75.4-z 



In modern metallurgical parlance, a martensitic trans- 
formation has come to denote any diflFusionless, struc- 
tural transformation resulting in a long-lived metastable 
phase with a high degree of short-range order. For in- 
stance, when adiabatically cooled below Tc = 910°C, 
Fe undergoes a first-order structural transition from an 
FCC solid (austenite) to an equilibrium BCC solid (/er- 
rite). A faster quench produces instead, a rapidly trans- 
formed metastable phase called the martensite. On nu- 
cleation, martensite 'plates' (~ 1/im) grow with a con- 
stant front velocity {r^lO^cm s~^) ^ ; fast compared to 
typical atomic diffusion speeds. The plates consist of an 
alternating array of twin BCC crystals along the equa- 
torialplane. Despite the considerable amount of empir- 
ical Q and theoretical work a unified theoretical 
approach addressing both kinetics and morphology, and 
capable of describing the metastable martensite as well 
as the equilibrium ferrite, has not emerged. In this Let- 
ter, we present a mode-coupling theory for the nucleation 
and growth of a product crystalline droplet within a par- 
ent crystal. We show that for slow quenches, the droplet 
grows diffusively as an equilibrium ferrite inclusion, while 
for fast quenches, the droplet grows ballistically, as a 
martensite having twinned internal substructure, with a 
speed comparable to the sound velocity. 

Our mode-coupling dynamics involves "slow" modes 
that change over time scales corresponding to the prop- 
agation of the nucleation front. The slow modes of a 
solid undergoing a first-order structural transformation 
are the displacement field u(r, i), and the momentum 
density g(r,t). Imagine, however, a droplet of the prod- 
uct nucleating within a parent crystal at time t = 0. It is 
clear (Fig. 1(a)), that an atomic mismatch is generated at 
the parent-product interface as soon as the nucleus is 
formed. This mismatch appears as a discontinuity in the 



normal component of the displacement field across the 
parent-product interface |^ Au • n = 0, and leads to a 
compression or dilation of the local atomic environment. 
Fig. 1(b). Since the martensite front velocity is large 
compared to atomic diffusion times, 0(r, t) appears as a 
new slow mode (vacancy field), which relaxes diffusively 
over a time scale r^. 




FIG. 1. (a) Inclusion of a product (rhombus) in a par- 
ent crystal, showing the discontinuity in u across the parent- 
product interface at t = 0. (b) A 3-d plot of the initial value 
of (f){x,y) for the inclusion shown in (a). The four corners in 
(a) and (b) are clearly marked. 

A Langevin description of the dynamics of the transfor- 
mation, requires a free-energy functional which describes 
all intermediate configurations between an austenite and 
a ferrite. Thus the usual elastic free-energy functional of 
a solid, Fei, has to be augmented by an interfacial term, 
describing the parent-product interface. As we see below, 
such droplet configurations do not affect equilibrium be- 
haviour, though their effect on dynamics is significant. 

Nucleation of a droplet of the product of size L and 
interfacial thickness ^, results in a strained crystal. The 
free-energy functional |^] of this strained crystal is given 
by, 
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c(|R-R' + u(R)-u(R')|) 



(1) 



where {R} represent the lattice vectors of the parent 
crystal, u(R) are the displacement fields and c(|r — r'|) 
is the direct correlation function [|j of the liquid at freez- 
ing. The above expression is exact at T = ; corrections 
are of the order of the r.m.s. fluctuations of the atoms 
about their perfect lattice positions, which are small in 
the solid phase 

Since the range of c is of order ^, the sum over {R, R'} 
in Eq. (|l|), can be split into two parts, with {R, R'} being 
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on the same or on either side of the interface. Expand- 
ing u(R') about R in the first part, leads to the usual 
bulk elastic free-energy functional Fei (9|. However, as 
argued above, the assumption of continuity breaks down 
for the second part, and so such an expansion cannot 
be carried out. In the limit ^/L <^ 1, the interface can 
be parametrized by a sharp surface T (r) = 0, with R 
and R' lying infinitesimally close to F = 0. Across this 
interface, the normal component of u is discontinuous ; 
^iTO|R-R'Ho(u(R) - u(R')) • n = (u+ - u_) • n EE 0, 
where n is the unit normal to the interface. The discon- 
tinuity (f>{r), which is simply ^ (/o(r)|r=o — p)/P: is the 
local vacancy field (p and p are the local and the aver- 
age densities respectively). In the continuum limit, this 
leads to the free-energy functional (to leading order in 
the discontinuity), 

^ = Fei + ^ J ^ [(u+-u_).nf C5(r) , (2) 

where fio is the unit cell volume of the parent. The coeffi- 
cient 7 = Qq^^ Y1' £,^9ndnc{r), where the prime denotes a 
sum across the interface, is clearly the surface compress- 
ibility modulus of vacancies, whose magnitude is of the 
order of the bulk elastic moduli ; it is however dependent 
on the local orientation of the parent-product interface. 
Since the strains = (diUj + djUi)/2 are continuous 
across the interface, we replace the delta function in Eq. 

by a regulator [1 — exp(— ^ i9„ey )^], so that the 
final regulated free-energy functional which incorporates 
all the slow modes in the problem is, 

^ = i^ei + ^ y rfr (a„ e,, f . (3) 

At the initial time (t = 0), the transformed region 
(product) is simply obtained as a geometrical deforma- 
tion of the parent. Fig. 1(a), which fixes the initial value 
of (j) (Fig. 1(b)). Having created this discontinuity at 
the parent-product interface, it will diffuse over a time r^. 
The Langevin equation describing the linearised dyamics 
of the slow modes ui and 4> can be written as, 

p.. 5T 5T 

+ v-V0 = -V2^, (5) 

0(p 

where is the time derivative of the strain tensor e^ . 
The inertial term Ui (propagation of sound waves) and 
solid viscosity tensor 1/^^; (damping of sound waves |^) 
are obtained from the coupling of u with the momentum 
density g 

We study the consequences of this dynamical theory, 
in the simple setting of a first-order structural trans- 
formation from a square (austenite) to a rhombic (fer- 
rite) crystal in two dimensions. Our results, which 



can be easily extended to the tetragonal to orthorhom- 
bic transition (essentially a 2-dim square to rectangular 
), are of relevance to structural transformations in al- 
loys like In-Pb, In-Tl, Mn-Fe. This transformation in- 
volves a she ar-|- volume deformation, and so the strain 
order parameter has only one nontrivial component 

63 = {Uxy + Uyx)/2. 

We construct a bulk elastic free energy with three min- 
ima — one corresponding to the undeformed square cell 
(63 = 0) and the other two corresponding to the two 
variants of the rhombic cell (63 = ±eo). The free-energy 
functional, in dimcnsionless variables, is 

^ = / ae3'-e3* + e3^ + (Ve3)'+7(05„e3)' (6) 

Jx,y 

The three minima of the homogeneous part of F at 63 = 
(austenite) and 63 = ±eo = ±[(1 -h VI - 3a)/3]i/2 (fgj._ 
rite), are obtained in the parameter range < a < 1/3. 
The parameter a is the degree of undercooling T — Tc- 
The surface compressibility 7 for an interface between 
the square (63 = 0) and the rhombus (63 = eo) posi- 
tioned along the y-axis is given by 7 = (c"(Z) -I- c"{l + 
ef)l)) + c"{\l — eol\))/2, where c"{r) is the second der- 
itive of a typical direct correlation function for a two- 
dimensional fluid whose range is taken to be of the order 
of the distance between next-nearest neighbours of the 
parent square lattice of spacing I. 

There are two time scales relevant to our kinetics — 
the quench rate = {dT/dt)/T and the vacancy re- 
laxation time T0. Accordingly two extreme dynamical 
limits suggest themselves. When t ^ t^, the (p fields re- 
lax instantaneously (fast mode) to (f> — 0, its equilibrium 
value. The only surviving slow modes are Ui, which obey 
Eq. (HI), with the free-energy functional J- — F^i. Since 
we are interested in the growth of the product nucleus in 
the parent matrix, it is appropriate to rewrite Eq. (^ in 
terms of the broken symmetry mode 63 (r, t), which takes 
the form e3(r — R(t)), when the interface is sharp (R(t) 
is the position of the interface). The equation for 63 is 
purely dissipative, 63 = —SFei/Se^. Minimising F^i with 
respect to 63, subject to boundary conditions 63 = at 
|r| 00 and 63 = ±eo at |r| = 0, obtains a ferrite nu- 
cleus, growing as i? ~ i^/^ at late times ||l0| . 

On the other hand, when r ^ t^, the (j) fields are frozen 
in the frame of reference of the nucleating front. In this 
limit, (j) is nonzero only at the parent-product interface, 
and so v in Eq. (|^) can be interpreted as the local front 
velocity. In the frame front, (j) obeys a diffusion equation 
with the diffusion coefficient I?^ = a/t/t^. Recasting 
Eq. (H), in terms of 63, leads again to a purely dissi- 
pative equation with the free-energy functional J- given 
by Eq. (|6|). To determine the structure of the product 
nucleus, we use a variational ansastz for 63, consistent 
with boundary conditions mentioned above. Consider a 
rectangular nucleus of length L (along x) and width W 
(along y) divided into N twins (Fig. 2 (inset)). The TV — 1 
twin interfaces, all of thickness rj, can be parametrized 
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by ea (x), which for the i-th interface takes values — Cq 
at {i - l)L/N + 7^/2 < X < iL/N - ry/2 and eo at 
iL/N + r]/2 < X < {i + l)L/N - 77/2, connected by a 
Unear interpolation. The strain at the austenite-ferrite 
interface varies linearly between at a; < — and ±eo 
at X > 1^/2. The free energy E{L,W) for the nucleus 
reads, 

E{L, W) = [ AFLW + {21 ~ AF)m + WC) + 

(TV -!)(/- AF)Wr] + 2eo\L/^ + W/C) + 
W{N - l)Wh ] + {W/0[ 2L^N-^/i + 
L'^'niN'^ - N-'^) + Lif {1/2- -iN-^l A) - 
LQN-^LN-^ - C/4 - 77/2) ] + 
[2W^/2,-W'^S. + W(^/2] . (7) 

In Eq. (0), AF = ae^^ — cq"^ + gq^ is the difference be- 
tween the hulk free energies of the austenite and the fer- 
rite, while I = aeo^/3 — + eo^/7. The structure of 

the growing nucleus is obtained by minimising E{L, W) 
with respect to the order parameter profile 63 (both the 
amplitude and phase). Within our variational scheme, 
this amounts to minimising E{L, W) with respect to the 
interfacial widths 77, ^ and Q and the phase of the order 
parameter, N . Minimisation yields the structural rela- 
tion (Fig. 2), 



(8) 



The exponent, a ~ 1/2, with tiny deviations for large W . 
The exponent a is empirically known to lie between 
0.4 and 0.5. Note that our theory suggests that Eq. (||) 
holds at all times during the growth of a martensite as 
a consequence of local equilibrium. Our prediction for 
L/N can be verified from in situ TEM studies of growing 
martensite fronts. 

Twinning is a consequence of the (j) term which is con- 
fined to the parent-product interface. When t <^ t^, 
the quench nucleates an inclusion which initially grows 
as a ferrite. Further growth as a single-domain ferrite is 
discouraged since the energy cost at the interface is pro- 
portional to the square of the discontinuity in u. The 
growing nucleus, gets around this by creating a twin. 
Though this costs interfacial energy due to (Ves)^ (which 
is small), the contribution from is identically zero at 
the twin interface, due to symmetry. 

From the contour plots of the optimised free energy 
plotted in Figs. 3(a) and 3(b), one observes a minimum 
in W for large 7 ~ 0(1) (Fig. 3(b)). Since we find that, 
Wmin < L, the free energy is minimised by thin rect- 
angular strips reminiscent of 'hard' acicular martensitcs 
seen in Fe-Ni or Fe-C systems. For small 7 ^ 0(10^"^), 
there is no such minimum in VF. As we shall see below, 
growth along L is much faster than along W, and so the 
martensite traverses the entire extent of the sample along 
L, whereupon growth proceeds along W. Such single in- 
terface growth is indeed seen in 'soft' solids, like In-Tl 
and Au-Cd alloys [fll. 
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Fig. 2 Plot of L/N vs W^^^ for a twinned rectangular 
strip of length L and width W having TV twins, obtained by 
minimising E{L,W), for a — .01, 7 — .001. The inset shows 
the variational profile for ea, with interfacial widths ^ , ( and 
77. The -I- and — denote regions where 63 = ±eo respectively. 

The critical nucleus is obtained by locating the saddle 
point in each of these free ener gy surfa ces. Taking the 
curvature at the corners k = ^JtJ^jLW as a measure of 
the sharpness of the growing nucleus, we find the criti- 
cal nucleus to be diffuse for both 7 = 1 (k = 0.4) and 
7 = 10"'^ (k = 1.43). In situ TEM measurements on 
critical martensite nuclei, reveal such diffuse structures 
|]l2| . Once past the critical size, the martensite droplet 
grows. As growth proceeds the interfacial widths shrink 
and K rapidly reduces to a microscopic value. A well de- 
fined front velocity now emerges. We can determine the 
late time growth behaviour of the martensite droplet, by 
requiring that the rate of energy change, computed from 
the time derivative of E{L,W,'N ^ L/W), equals the 
energy dissipation from the local evolution of the order 
parameter 63 ||l^. Since 63 follows a purely dissipative 
dynamics. 



E {L, W,N L/W) ^ ly eldxdy 



(9) 



The time derivative of 63, can be easily computed from 
our ansatz, and leads to the following form for the en- 
ergy dissipation, I'eo'^ {L'^W/C + W^L/^). It is easy 
to determine the asymptotic solutions for L and W, 
from the resulting first-order, nonlinear ODE. Assum- 
ing L{t) = VLt" and W{t) = vwf^ as t ^ 00 (a,/3 > 0), 
and equating the dominant singular contributions on ei- 
ther side of Eq. (|), we obtain a = 1 and (3 = 1/2. The 
nonlinear ODE also gives the velocity ~ \fljv. Since 
7 is of the order of typical elastic moduli, the marten- 
site grows with a constant velocity close to the velocity 
of sound, in the direction perpendicular to the twinning 
plane. Parallel to the twinning planes, the front moves 
diffusively Q. 
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Before we close, we would like to point out some impor- 
tant features not included in our theory. When the sur- 
face compressibility 7 is large, the product might prefer 
(b) to generate dislocations at the parent-product interface, 




producing instead internal slip bands [gj. Moreover for 
solids with low thermal conductivity, transport of heat 
across the interface might significantly alter the shape of 
grains. We are currently working on these refinements to 
the theory. 

We thank S. G. Mishra and T. V. Ramakrishnan for 
discussions and encouragement. 
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Fig. 3 Contour plots of the minimized E{L, W) in the 
L-W plane for a = .01 and (a) 7 = 0.001 and (b) 7 = 1. 
A comparison of our value for rj with typical experimental 
values, yields a unit of length of 1 — 2 A. 

Thus the martensite nucleus grows in L with veloc- 
ity vl, till it collides with other growing nuclei, where- 
upon the high elastic energy barriers prevent coalescence 
psf . Growth in W proceeds slowly, upto a point where 
W = Wrain (for large 7). When 7 is small, growth in W 
proceeds unimpeded, unless pinned by other nuclei or im- 
purities. The solid has now got trapped as a metastable 
martensite and requires large thermal activation to trans- 
form to the equilibrium ferrite. The martensite can there- 
after reduce interfacial energy by deforming the crystal 
such that the discontinuity in u (and hence (/)) vanishes. 
This will, however, generate long-range stress fields 
In addition, growth can terminate, if is made to relax 
(say, by raising T), leading to thermal arrest of marten- 
sites 0. 

A more detailed account will explore a 'morphology 
phase-diagram' as a function of kinetic and structural 
parameters. In this paper, we have confined ourselves to 
a simple variational ansatz for the shape. Indeed, hard 
martensites like Fe-Ni are 'lens shaped' plates. We see 
however, that the interfacial energy at either end of our 
rectangular strip can be further reduced by decreasing 
W{x = L) and W{x — 0), thus forming a 'lens'. In 
three dimensions, the possibility of several twin variants 
would lead to interesting stacking patterns. Moreover, 
the fascinating phenomenon of shape memory can be 
studied by coupling to an external stress field. 

Conventional analysis of the kinetics of martensites 
minimises the elastic energy, subject to the con- 
dition that u is continuous across the parent-product in- 
terface. This ad-hoc boundary constraint gives rise to 
long-range stress fields. The analysis of the dynamics is 
complicated by the fact that the boundary constraints 
are moving, allowing for analytic solutions only in 1-dim 
1^. Our approach does not suffer from this handicap, 
since the (j) field, which lives at the interface, naturally 
leads to the boundary constraint at late times. 
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